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Abstract. Using a two-channel model, we show that the number of closed- channel molecules in a two- 
component Fermi gas close to a Feshbach resonance is directly related to the derivative of the energy 
of the gas with respect to the inverse scattering length. We extract this quantity from the fixed-node 
Monte Carlo equation of state and we compare to the number of closed-channel molecules measured in 
the Rice experiment with lithium [Partridge et al., Phys. Rev. Lett. 95, 020404 (2005)]. We also discuss 
the effect of a difference between the trapping potentials seen by a closed-channel molecule and by an 
open-channel pair of atoms in terms of an effective position-dependent scattering length. 

PACS. 03.75.Ss Degenerate Fermi gases - 67.90.+Z Other topics in quantum fluids and solids 



1 Introduction 



It is now experimentally possible to prepare two-component 
atomic Fermi gases at low temperature with a fully con- 
trolable interaction strength, in a regime where the inter- 
action range is negligible and the atomic interactions are 
characterized by a single parameter, the s-wave scattering 
length a between two opposite "spin" atoms. The value 
of a can be adjusted at will thanks to a magnetically in- 
duced Feshbach resonance. The weakly interacting limits 
kpa — > 0^ and kpa — + 0^, where kp is the Fermi wavevec- 
tor of the gas, correspond respectively to the BEC limit (a 
Bose-Einstein condensate of dimers, observed in [li2,3,!4, 
[5]) and the BCS limit (a "condensate" of Cooper pairs ap- 
proximately described by the BCS theory). Furthermore, 
the experiments can access the so-called crossover regime 
between BEC and BCS, where the gas is strongly inter- 
acting (fcF|a| ^ 1), which includes the celebrated unitary 
limit A;_F|a| = oo, where the gas acquires fully universal 

In this context, early theoretical studies [l8lfT9ll20] put 
forward a many-body Hamiltonian which accurately mod- 
els the microscopic two-body physics of the Feshbach reso- 
nance. The interaction potential is described by two chan- 
nels, an open channel and a closed channel. The atoms ex- 
ist in the form of fermionic particles in the open channel 
and in the form of bosonic short range molecules in the 
closed channel. Two atoms may be converted into a short 
range closed-channel molecule and vice versa due to the 



interchannel coupling. The corresponding Hamiltonian is 
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V'fc(r) 
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ri+r2 r3+r4\ t(,^W(r2)^^(r3)^^(r4) (i) 



where the atomic fields ipa (r) obey the usual fermionic an- 
ticommutation relations and the field 4'b (r) describing the 
closed-channel molecules obeys the usual bosonic commu- 
tation relations. A gives the coupling between the closed 
and open channels, with a short range cut-off function 
x(r), of range &, and such that f d'^rxM = 1. A typi- 
cal value for b is in the nanometer range Q|. The coupling 
constant go represents the background atomic interaction 
in the open channel, which is conveniently modeled by a 
separable potential with the same cut-off function x- The 
energy Ei, is the energy of a closed-channel molecule, in 
the absence of the coupling A, counted with respect to 
the dissociation limit of the open channel. This energy 
is adjusted with a magnetic field B, using the fact that 
the different magnetic moments in the open and closed 
channels lead to a differential Zeeman shift. The resulting 



^ 6 can be estimated by the van der Waals length, that is 
the length that one forms with fi, m and the Ce coefficient of 
the van der Waals interaction 1211. 
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effective magnetic moment is 



fJ-b 



dB ■ 



(2) 



C/(r) and Ub{r) are the trapping potentials experienced by 
the atoms and the closed-channel molecules respectively. 
On the contrary, recent Monte Carlo simulations of 
the many-body problem in the BEC-BCS crossover sim- 
ply use single channel model Hamiltonians [22,23,24,25i, 
[26 1 127 ( 128] . In most unbiased Quantum Monte Carlo simu- 
lations, a lattice model is used: 



^i = E^' E^iw 



. <7=Ta 



ft2 

- — A + U{r) 
Zm 



Mr) 



-goV'|(r)V'I(r)^x(r)^T(r)}- (3) 



Here b is the lattice spacing (which coincides with the 
"range" of the discrete delta interaction potential), A is 
a discrete representation of the Laplacian, and go is a 
coupling constant adjusted to have the desired scattering 
length: As was shown e.g. in [24l[25l[29ll30] . 



1 
9 
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— Tr/b.Tr/b] 



d^k 1 

3 (2^2^ 



(4) 



where g — Awfi^a/Tn is the effective coupling constant and 
eic the energy of a single particle with wavevector k on the 
lattice. In fixed-node Monte Carlo calculations, one rather 
uses an interaction potential in continuous real space, e.g. 
a square well potential [3T 1 l32 l [33 l l34 p 35 p 36j . 

There is now a consensus about the fact that, in the 
zero-range limit, there is a convergence of predictions for 
the two models Hi and H2 , in particular for the equation 
of state of the gas. The zero-range limit means, for single 
channel models, that kpb <C 1, & <C A, and b ^ \a\, where 
kp is the Fermi wavevector and A is the thermal de Broglie 
wavelengthjj For two-channel models, one has to impose 
the same conditions, not only for the interaction range b 
but also for the effective range rg (one indeed has \re\ 3> b 
on narrow Feshbach resonances [371138] ): we shall see in 
appendix 1X1 that one also needs to impose a condition on 
the background scattering length, and that all conditions 
are satisfied in the experiment [5]lfl 

One can however consider observables which exist in 
reality and are included in the two-channel model, but 
which are absent in a single channel model. The most nat- 
ural example is the number of closed-channel molecules, 



^ The assumption b <^ \a\ is not necessary on the BCS side 
if one restricts to macroscopic observables such as the total 
energy, but it is necessary for microscopic observables such as 
the ones considered in section |3TT] 

^ A simple formal definition of the zero-range limit for the 
homogeneous gas is to take the limit of vanishing density kp — > 
while tuning the magnetic field B in such a way that iKkpo) 
remains constant, and keeping constant all parameters of the 
Hamiltonian other than B ; one then also has to take the limit 
T ^ so that T/Tf remains constant. 



which was recently measured at Rice using laser molecular 
excitation techniques [S]. 

In this paper, we show that the equilibrium number 
of closed-channel molecules Nb is directly related to the 
equation of state of the gas, more precisely to the deriva- 
tive of the gas energy (or free energy at non-zero tem- 
perature) with respect to 1/a. Paradoxically, we can thus 
use the equation of state calculated in 32J within a sin- 
gle channel model, in order to predict N^ in the crossover 
and compare to Rice measurements, see section [2l Since 
the derivative of the energy with respect to 1/a is related 
to other observables involving atomic properties only [39l 
140] , it is also possible in principle to access N^ by a pure 
atomic measurement, see section [3j We conclude in sec- 
tion gj 



2 Prediction for the number of 
closed-channel molecules 

2.1 General result 

We first suppose that the system is in an arbitrary eigen- 
state |'0) of the two-channel Hamiltonian _ff2, with eigen- 
ergy E. From Hellmann-Feynman theorem we have 



dE dH2 



(5) 



The only magnetic field dependent part of the Hamilto- 
nian H2 is the bare closed-channel molecule energy Eb{B) 

so that 

dE 

dB - ^^""^ (^) 



where 



iV,= d\{i^l{r)Mr)) 



(7) 



is the mean number of closed-channel molecules. 

We then eliminate the magnetic field by using as a pa- 
rameter the scattering length a experienced by the trapped 
atoms rather than B. In the general case where the trap- 
ping potential C/fc(r) for the closed-channel molecule differs 
from twice the trapping potential U{r) experienced by an 
open-channel atom, the scattering length of two atoms 
around point r depends on r, see appendix [BJ To lift the 
ambiguity we thus take for a the scattering length in the 
trap center taken as the origin of coordinates, r = 0. We 
then rewrite ([Sj) as 



Nh 



dE d{l/a) 1 
d{l/a) dB ^b 



(8) 



Assuming that fib is magnetic field independent, a as a 
function of B can be calculated explicitly for H2 by solving 
the zero-energy two-body scattering problem in free space 
(see appendix [A]) . and by including the energy shifts due 
to the trapping potentials Ub{0) and U{0), as explained 
in the appendix [BJ The function takes the form 



ai,iB-SBo) 



(9) 
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where we have introduced the function giving the free 
space scattering length as a function of B, 



afe(B) = flbg 1 - 



AB 
B-Br 



(10) 



Here Obg is the background scattering length, Bq is the 
location of the Feshbach resonance in the absence of trap- 
ping potential, SBq = [2C/(0) - Ub{0)]/fib the shift of the 
resonance location due to the trapping potential, and AB 
the resonance width. 

Equation ([8]) is directly applicable at zero tempera- 
ture. At non-zero temperature, one has to take a thermal 
average of ([5]), keeping in mind that the mean value of a 
derivative is not necessarily equal to the derivative of the 
mean value. In the canonical ensemble one can check the 
exact relation WU 



dE 

d{l/a) 



dF 
d{l/a) 



d{E) 
d{l/a) 



(11) 



where (. . .) is the thermal average, the derivative of the 
free energy F is taken for a fixed temperature T, and the 
derivative of (E) is taken for a fixed entropy S. We thus 
have in the canonical ensemble 



(^;>> = 



f d{E) 
\d{l/a) 



d{l/a) J_ 

^J'b 



dB 



(12) 



2.2 Analytical results in a trap in limiting cases 

We now restrict to a spin balanced gas where the number 
of particles is equal to N/2 in each spin component. We 
also assume that the gas is at zero temperature (or at 
temperatures much smaller than the Fermi temperature 
Tp), in a harmonic trap 



U{v) = C/(0) + im 



E 



2„2 



w^n 



(13) 



in the macroscopic regime where ksTp is much larger than 
the oscillation quanta fkUa. The usual definition of the 
Fermi temperature Tp and of the Fermi wavevector kp is 
then 

ksTp ^ ^ ^ (37V)i/3to (14) 

Zm 

where iu is the geometric mean of the three oscillation 
frequencies Wq. In this macroscopic regime, we rely on the 
local density approximation to calculate dE/d(l/a) for the 
trapped gas: This is expected to be exact in the large N 
limit, and to already hold for the Rice experiment [S] since 
ksTp/ihuj^) > 9. 

The key ingredient of the local density approximation 
is the equation of state of the interacting homogeneous 
gas. At zero temperature, and in the zero-range limit de- 
tailed in the introduction, the equation of state is expected 
to be a universal function of the gas density and of the 
scattering length, independent of the microscopic details 
of the interaction. Furthermore, the spatial dependence of 



the scattering length in presence of the trapping poten- 
tial, due to Ub{r) ^ 2 C/(r), is expected to be negligible in 
the zero-range limit, as argued in the appendix iBl Under 
these assumptions, the local-density approximation for the 
trapped gas leads to the universal form 



dE 
d{l/a) 



NT 



kpa 



(15) 



where .F is a dimensionless function, expressed in the ap- 
pendix [C] in terms of the dimensionless universal function 
/ giving the energy per particle e'^°"^ of the homogeneous 
interacting gas: 



^hom 



^lioin , 



1 






(16) 



Here eo°™ and k^™ are the energy per particle and the 
Fermi wavevector of the ideal Fermi gas with the same 
total density n as the interacting gas: 



7,hom 

h.p 



horn 
^0 



1/3 



(BTT^n) 

3 h^ {k^°"'y 
10 m 



(17) 
(18) 



In the EEC regime < kpa < 1 we have [1^1431144] 

1 \ 5 



uhom , 



3(fc^°ma)2 



(19) 



5 
IStt 



fc^°"^c 



:(4™^arf)3/2 



where the first term comes from the dimer binding energy 
}t? /(mo?), the second term is the mean field interaction 
energy among the dimers, proportional to the dimer-dimer 
scattering length a^, and the last term is the (bosonic) 
Lee- Huang- Yang correction. The solution of the 4-fermion 
problem gives [45u46j 



ad — 0.60a. 
The local-density approximation then leads to 
1 \ _ 1 5^/5 ^arf^2/5 



T 



kpa 



kpa 212/5 . 7 V 



{kpa) 



7/5, 



(20) 



(21) 



where we have omitted the Lee-Huang- Yang correction. 
In the BCS hmit < -/cfo < 1, we have PUiipTl 



uhom , 



1 



r,hom 

4(11 -21n2) 
21^;^ 



(A:|°'"a)^ 



(22) 



where the first term is the ideal gas result, the second 
term is the Hartree-Fock mean field term, and the last 
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term is the (fermionic) Lee-Huang- Yang correction. This 
leads toQ 



T 



kpa 



512 

' 9457r2 
/ 256 
1^35^ 



[kpaf 



1+ 



63 -f 189 In 2 
1024 



kpa 



(23) 



Around the unitary limit \/{kpa) ^ wc use the ex- 
pansion 



/ 



Uhom 
■up L 



^ 7,ho 



which gives for the trapped gas: 

2^ C 



^(0) 



52-7-7r^i/4 
where we took the estimates 



0.28 



0.41 
0.95 



(24) 



(25) 



(26) 
(27) 



The estimate for ^ is obtained by averaging the predic- 
tions of [32] (^ = 0.42(1)), of [26| (^ = 0.449(9)), of [33] 
(C = 0.42(1)) and of |^ (^ = 0.37(5)). The estimate for 
C, is obtained from the calculation of the short range pair 
correlation of [34] , see subsection 13.11 and is close to the 
value C ~ 1.0 that may be extracted directly from the 
values of / in f32l by approximating the derivative by a 
two-point formula. 13 



2.3 Comparison to Rice experiment 

At Rice [5 the quantity Z = Nb/{N/2) was measured for 
a lithium gas, by resonantly exciting the closed-channel 
molecules with a laser that transfers them to another, 
short-lived molecular state; the molecule depletion reflects 
into a reduction of the atom number that can be measured. 
We now compare the Rice results to our prediction ([8|) . 
We first insert (|9I10|) into ([8]) which gives 



Nb = NkpR^T 



1 
kpa 



1 



obgy 

a J 



(28) 



* For a small positive scattering length < kpa <^ 1, 
Eas. p2l23p also apply to the atomic gas state; since the in- 
teraction potentials considered here have a two-body bound 
state for a > 0, this atomic gas state is only metastable [49], 
the ground state being the BEC of dimers considered previ- 
ously [cf. Eq.([2l|]. 

^ The values (|26l27p of ^ and ^ can be inserted into the ex- 
pressions [SOllST] for the derivatives of the frequencies of the 
hydrodynamic collective modes with respect to l/(fc_Fa) taken 
at unitarity. The resulting slope is compatible with the exper- 
imental data for the radial mode [121152] . while the agreement 
with the data for the axial mode [53] is only marginal. The 
data from [12] agree very well for all positive values of a with 
the theoretical result [54) obtained from the fixed-node Monte 
Carlo equation of state [32], while in [521153] finite-temperature 
effects may play a role [12II54II55] . 
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Fig. 1. Derivative —dE/d{l/a) of the energy _B of a trapped 
two-component Fermi gas, where a is the s-wave scattering 
length. Solid line: Theoretical prediction combining the Lee- 
Huang- Yang formulas (]19l22p in the weakly interacting regimes 
with an interpolation of the fixed-node Monte Carlo results of 
[321134] in the strongly interacting regime (see appendix JD)) . 
Dashed lines: Analytical predictions in the weakly interacting 
regimes, on the BEC side ([2T} and on the BCS side (f23)l . Cross: 
Result ([25p at unitarity. Circles (with error bars): Experimen- 
tal measurement of the number Nb of closed-channel molecules 
in a lithium gas at Rice [5], combined with the present theory 
linking dE/d(l/a) to Nb, see ([28p . (a): Linear scale on the ver- 
tical axis, (b): Logarithmic scale. —dE/d{l/a) is expressed in 
units of Nh^kp/m. All the theoretical predictions are obtained 
in the local-density approximation. 



where the function T is defined in ([T5|) . The length i?* 
depends on the interchannel coupling and is related to 
the width of the Feshbach resonance as: 



i?* 



mahgl^b^B 



(29) 



For the Bq — 834G Feshbach resonance in ^Li, we have 
fib ~ 2/iB where fiB is the Bohr magneton J21], AB = 
-300G, abg = -1405ao where gq = 0.0529177 nm is the 
Bohr radius [53]; this gives i?* — 0.0269 nm. From the 
Rice data for Nb/N, using ([28| . we thus get values of T, 
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i.e. of dE/d{l/a), that we compare in Fig[T] to theoreti- 
cal predictions, given in the local-density approximation 
either by the expansions (|21l23p in the weakly interacting 
regime, or by ([25|) at the unitary limit, or by an interpola- 
tion of the Monte Carlo results of 532','34] in the strongly 
interacting regime (see appendix [C] for details) . The 
agreement is satisfactory, the experimental points being 
scattered around the theoretical prediction. In particular, 
a possible shift SBq of the resonance location due to the 
fact that C/fc(r) 7^ 2U{r), see (j9|), an effect not explicitly 
included in [5, for the calculation of l/(kpa) as a function 
of B, and which would manifest itself in horizontal shifts 
of the experimental data with respect to theory, does not 
show up in Fig[T] For the considered broad resonance of 
lithium, this is not surprising, see appendix [Bl 

Our approach even allows to predict the time evolu- 
tion of the particle number, in presence of the depleting 
laser, under the assumption that the molecular depletion 
rate is so weak that the equilibrium relation (|28p holds 
at all times. Neglecting the slow off-resonant free-bound 
photoassociation process [5] , we have 






(30) 



where >!7^/7 is the effective decay rate of the closed-channel 
molecule, the Rabi frequency fl and the spontaneous emis- 
sion rate 7 being defined in [5]I3 

Simple expressions can be obtained in limiting cases: 
In the BEC limit, the dominant term of the expansion 
([2T|1 gives 

N{t)^N{0)e-~^°*; (31) 



in the unitary limit, (j25p gives 



Nit) 



N{0) 



(i + rot/er^ 

and in the BCS limit the dominant term of (P5|) gives 



(32) 



N{t) 



N{0) 



{i + rot/2y 



(33) 



Here, /q = —[dN{0)/dt]/N{Q) is the inital atom loss rate, 
which can be expressed in terms of the initial atom number 
7V(0) using Eqs. (|28I30|) . We see that in general N{t) is 



^ The error bars in Fig. 1 include the systematic theoretical 
uncertainty resulting from the leading order correction to (|10|) 
given in [56], but are largely dominated by the experimental 
error bars from [5] because we did not include the three ex- 
perimental points from [^ which are deep in the BEC regime 
[l/{kFa) > 5]; in this regime, two-body theories which go be- 
yond our two-channel model give good agreement with the Rice 
measurements [5l l21|[57] . 

''' Eq. (|30p can be obtained from the following master 
equation for the density matrix p of the gas: dp/dt = 

{ny'y)Jd'r [Mr)p^l{r) - U4MMr^),P}] + [H2,p]/{ih). 
N in (|30|l stands for TVf + 2Nb, where A^'f is the number of 
fermions; thus A'^ ~ TVf in the regime considered in this paper. 



not an exponential function of time. For the last point on 
the BCS side in [5], the (unpublished) data for N{t) are 
better fitted by Ea. (|33p than by an exponential; however 
the value of Fq resulting from this new fitting procedure 
does not differ significantly from the published value [5Sj . 



3 Related observables 

In this section, contrarily to subsections 2.2 and 2.3, we 
return to the general case and we do not assume that 
the temperature is ^ Tp and that the numbers of atoms 
in each spin component are equal. For concreteness we 
take the lattice model ^, even if the reasonings may be 
generalized to a continuous space model and to a two- 
channel model. 



3.1 Short range pair correlations 



As was shown in [391140] and recently rederived in [59] , in 
the zero interaction range limit, the derivative of the gas 
energy with respect to the inverse scattering length can 
be expressed in terms of the short range behavior of the 
pair distribution function of opposite spin particles: 



d{E) 
d{l/a) 



47rr 



d'^Rlimr^gt|(R + r/2,R-r/2) 

r—>-0 

(34) 
where the pair distribution function is given by 

5Ti(ri,r2) = (V'|(ri)V|(r2)Vi(r2)V'T(ri))- (35) 

Since we are dealing here with a finite range interaction 
model, we have taken a zero-range limit in (|34| . introduc- 
ing the notation 

lim = lim lim . (36) 

Note that the order of the limits is important. The limit 
in (|34l) is reached for a distance r much smaller than 
|a|, the mean distance between particles, and the ther- 
mal de Broglie wavelength A, but still much larger than 
the interaction range b. 

We thus see that a measurement of the pair distri- 
bution function of the atoms gives access, using ([8]), to 
the number of closed-channel molecules. Experimentally, 
the pair distribution function was measured via the pho- 
toassociation rate in a ID Bose gas [60ll61[ . or simply by 
dropping the gas on a detector with high enough spatio- 
temporal resolution |62,63,64J . The possibility of measur- 
ing the pair distribution function in the BEC-BCS crossover 
was studied in [341165] . In the unitary limit, by inserting 
into Eq. fM]) the value of limr^g^i calculated in [31] with 
the fixed-node Monte Carlo technique, we get the value 
([27]) for the parameter ( defined in ([24]) . 

The relation ([M]) is the three-dimensional version of 
the relation obtained in [66] for a one-dimensional Bose 
gas with contact interaction, using the Hellmann-Feynman 
theorem. We now show that the Hellmann-Feynman tech- 
nique also provides a simple derivation of (|34p . Taking the 
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derivative of the gas eigenergy E with respect to the ef- 
fective coupling constant g — ATrfi^a/m, for a fixed lattice 
spacing b, the Hellmann-Feynman theorem gives 



dgo \r^ ,3 / N 



(38) 



where go is a function of g as given in ([4]), and where we 
used the relation (dE/dg) — {d{E)/dg)s [cp. pTjl ], In a 
gas, with an interaction range much smaller than the mean 
distance ^ 1/kp between particles and than the thermal 
de Broglie wavelength A, the pair distribution function is 
dominated by two-body physics if |ri — r2| is much smaller 
than 1/kp and A (and also much smaller than a in the 
BEC regime), a property contained in the Yvon '67' and 
Waldmann- Snider J38J ansatz of kinetic theory, which also 
appeared in the context of quantum gases [MUBMTnilTT] : 

gn(.ri,r,)^c(^^^±^^ |</'(ri-r2)P (39) 

where (j) is the two-body zero-energy free space scattering 
state for two atoms, here in the spin singlet state. With 
the lattice model we find El 



W")!^-^- 



so that Eg. ((55)) simplifies to 



m. = ?'=-<"- 



(40) 



(41) 



On the other hand, at interparticle distances much larger 
than b, the zero-energy scattering wavefunction (j) behaves 
as 1 — a/lri — r2| so that 



limr^5t|(R + r/2, R - r/2) = a^C(R). (42) 

r— >0 



Eq.dHI) thus leads to 



3.2 Tail of the momentum distribution 

As was shown in [?D] and recently rederived in [55], in 
the zero-range limit, the momentum distribution of the 
gas has a large momentum tail scaling as 1/k^, with a 



An elegant derivation is to enclose two atoms in a ficti- 
tious cubic box with periodic boundary conditions of volume 
L'^, with L 3> |a|, b. To leading order in a/L, the wavefunction 
is (pboxiri — r2) = (t>{ri — r2)/L^^^ with an energy given by the 
"mean-field" shift -Bbox = g/L"^ ■ The Hellmann-Feynman the- 



orem (|37p then gives the result (|40|) . An alternative derivation 
is to directly calculate (j) in Fourier space [see (|47p ] which gives 
0(r) = l-<?o0(O)/[_^/,^^/„3[d'k/(2^)^]exp(ik-r)/(26k); this, 
combined with (J4]), gives (|40p . 



coefficient common to the two spin states and proportional 
to -dE/d{l/a): 



lim k n„{k) 

k — *+oc 



47rm / d{E) 



ft2 \d{l/a) 



where 



lim = lim 

k — >+oo k — * + oc 



lim. 



(43) 



(44) 



a =1 or I and the momentum distribution is normalized 

as 

d-^k 

.^^n4k) = {N,), (45) 

where (N^) is the mean number of particle in the spin 
state a. When combined with ([5]), this relation reveals that 
the number of closed-channel molecules may be deduced 
from a measurement of the momentum distribution of the 
atoms. 

The momentum distribution in a two-component Fermi 
gas was actually already measured [72ll73l[74] . simply by 
switching off the interaction and measuring the cloud den- 
sity after a ballistic expansion. Two practical problems 
however arise, the effect of the non-zero switch-off time of 
the interaction [721175] and the signal to noise ratio in the 
far tails of the distribution; we have therefore not been 
able to extract the value of dE/d{l/a) from the published 
data for ndk). The value of dE/d{l/a) has not been ex- 
tracted either from the fixed-node Monte Carlo calcula- 
tions of Tig. (A;) performed in [76] . 

The relation ([^5]) is a three-dimensional version of the 
relation derived for a one-dimensional Bose gas in [77] . We 
now present a very simple rederivation of (|43|) . 

In a gas with short range interactions, i.e. b <C 1/kp 
and 6 ^ A (and 6 <C a on the BEC side) , the tail of the 
momentum distribution is dominated by binary collisions 
in the spin singlet channel between a particle with a large 
momentum k and a particle with a momentum ~ — k, so 
that at large k 

ri,(k)~S|^(k)|2 (46) 

where B is momentum-independent and (/)(k) is the Fourier 
transform J^r b^e~''^''^(f>{r) of the zero-energy two-body 
free space scattering state 0(r). For the contact interac- 
tion in the lattice model the two-particle Schrodinger's 
equation takes the form [50] 



2ek</'(k) -f 5o0(r = 0) = 0, 



(47) 



where ek is the free wave dispersion relation in the lattice 
model, so that (j>(k) drops as l/ck for large k, and the 
associated kinetic energy diverges for 6 ^ 0. Eq. fi^ im- 
plies that the kinetic energy of the gas is proportional to 
the kinetic energy of the zero-energy two-body scattering 
state in the zero-range limit: 



E^kinfeas] ~ 6-SkinM. 



(48) 



Since the total energies have a finite limit for 6^0, we 
can replace the kinetic energies by the opposite of the 
interaction energies in the above formula: 



Eint[gas]~S-Si„t[0]. 



(49) 
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(50) 



The interaction energy of the gas is go ( J2r ^^"^j V'lV'i. V"? ) i 
equal to {d{E) / dg) s go{dg / dgo) according to (P7|) . The in- 
teraction energy of </> is simply go 1 0(0) P, equal to go{dg/dgo) 
according to (PP]). The common factor goidg/dgo) simpli- 
fies and 

It remains to relate the coefficient B to limk~>+ook'^na-{k) , 
using (|46p . In the zero-range limit, the scattering state 
0(r) tends to 1 — a/r, so that, by Fourier transform, 0(k) 
tends to {2TT)^6{k) - Ana/k'^. We then get (|i5)) . 

4 Conclusion 

We have calculated the number Ni, of closed-channel mole- 
cules in the whole BEC-BCS crossover for a two-component 
Fermi gas. Our result is in satisfactory agreement with 
the experimental results from Rice [5]. The expression 
([8]) for Nb is proportional to the quantity —dE/d{l/a) 
where E is the energy of the gas; this quantity is univer- 
sal in the zero interaction range limit. At unitarity we find 
Nb/N — kpR^!F{Q) where .F(O) is a universal constant 
given in (|25p and the length i?* is related to the reso- 
nance parameters by (|29p : in the zero range and zero ef- 
fective range limit where kpb <^1 and fci?|re| ^ 1 we have 
kpR* < 1 [see Eq.([65l)], we thus confirm that Nb -^ N 
i.e. the gas mainly populates the open channel. In the 
BCS limit, —dE/d{l/a) and thus Nb are determined by 
the Hartree-Fock mean field energy, and do not tend ex- 
ponentially to zero with l/(fci?|a|), contrarily to the pre- 
diction of usual BCS theory [S^. This is also related to 
the fact that the BCS theory does not predict the correct 
short-distance pair correlations [78] . 

The quantity —dE/d{l/a) is related to the short-distance 
behavior of the pair correlation function and to the large-fc 
tail of the momentum distribution, as discovered by Tan 
[39ll40] and rederived in [59] and in this paper. The quan- 
tity —dE/d{l/a) is also proportional to the average ra- 
diofrequency shift [7^150] . The general idea that various 
short-range two-body quantities such as the number of 
closed-channel molecules, the interaction energy of the gas 
and the average radiofrequency shift are the product of the 
value of the considered observable in the zero-energy two- 
body scattering state and a single universal many-body 
quantity was given by Leggett [71] . Our work shows that 
the Rice experiment [5] is the first direct measurement of 
this fundamental quantity in the BEC-BCS crossover. 

It would be interesting to study the effect of a fi- 
nite temperature. Experimentally, this could be achieved 
by measuring the number of closed-channel molecules to- 
gether with the density profile of the trapped gas, the 
released energy ,6,4 , or the entropy [I4jll5j. Theoreti- 
cally, Eq. fT^ can be used when more Monte Carlo re- 
sults at finite temperature will become available. A use- 
ful tool in this context is the virial theorem: at unitarity 
it directly gives the total energy in terms of the density 
profile pT|fT5,81,82 while for a finite a it also involves 
[d{E)/d{l/a)]s [83,i84j. 
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A Two-body analysis of the two-channel 
model in free space 

A.l Scattering amplitude 

In this appendix we briefly discuss the two-body scat- 
tering properties for our model Hamiltonian H2 in free 
space, that is in the absence of trapping potential {Ub{r) = 
0, C/(r) = 0). The scattering problem may be considered 
in the center of mass frame so that we take the two-body 
scattering state 



/34|0) -t- 



_fk_ 



A{k)al^al^\0), (51) 



with the unknown amplitude f3 in the closed channel and 
the unknown momentum dependent amplitude ^(k) in 
the open channel. The annihilation operators for the bosonic 
molecules, 6k, and for the fermionic atoms ak^, cr =t or |, 
are the Fourier transforms of the corresponding field op- 
erators ipb{j^) and ipa{^), where the Fourier transform of a 
function /(r) is /(k) = J d^rexp{—i'k ■ r)/(r). Inserting 
the state ([ST]) in Schrodinger's equation (if 2 — ^)|V') = 
and projecting respectively on the molecular subspace and 
the atomic subspace leads to 



{Eb ~E)f3 + A-f^0 (52) 
^E]A{k)+x{k){(3A + goj)^0 (53) 



h^e 



where we used the fact that the function x(r) is real and 
even, and where we introduced the amplitude 



1 



_fk_ 
{2nf 



x(k)A(k). 



(54) 



We now specialize to a state corresponding to the scatter- 
ing of a spin ] atom with incoming wavevector kg onto 
a spin [ atom with incoming wavevector — ko, so that 
E = lT?ko/m is the energy of the incoming state. Elim- 
inating (3 in terms of 7 with ([52]) and using ([53]) leads 
to 



A(k) = (27r)3,5(k-ko) 



-lx{k)[go + A^/{E~Eb)] 
E + iO+- h?k'^/m 



(55) 



where we have singled out the incoming state, and the 
scattered component is guaranteed to correspond to a pure 
outgoing wave in real space by the introduction of +iO^ in 
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the denominator. Injecting ([55|l in ([54|) gives a closed lin- 
ear equation for 7. From the known large r expression for 
the kinetic energy Green's function we extract the scat- 
tering amplitude 



fk 



-mx(ko)V(4^S2 



[50 






-1 _ r d-'>k 



X"(k) 



(56) 



We now consider the effective range re defined by the 
low-g expansion —l/fiq = a~^ ^ q + \<f^e + 0{q^). Its 
general expression for our model is 



26^ 
a 



(64) 



Of particular interest is the value of the effective range 
right on resonance, 



Since we are considering a s-wave Feshbach resonance, the 
function x(r) is rotationally invariant, a property that we 
have used in ([SS]) . 

We now perform explicit calculations of various quan- 
tities characterizing the low energy scattering, such as the 
scattering length and the effective range. To this end, it 
is convenient as in [38j to perform the specific choice 



46 

where the length i?* is always non-negative, 

2 



i?* — 



A 



27rfe^o 



(65) 



(66) 



X(k) = exp(-fc^672) 



(57) 



and to extend ((56|) to negative energies, setting fco = 
iq,q > and E = —fi^q^/m. One then obtains an explicit 
expression 

1 

Jiq 

(58) 

where the last term —q is imposed by the optical theorem 
and erf is the error function. 

The inverse of the scattering length is obtained by tak- 
ing g ^ in ((58)) : 



p-9^b= 


ATTfi^/m ^ 1 


+ gerf (g6) - 








q 




190 Ei+K^qi/m "^'N 





This definition of i?* agrees with Eq. 

A very convenient rewriting of the scattering ampli- 
tude fiq is obtained by taking as independent parame- 
ters, in addition to the van der Waals length 6, the two 
lengths UhgjR* rather than the bare parameters go, A of 
the Hamiltonian. The quantity Ej, — A'^ /go appearing in 
(|58)l is conveniently rewritten using the quantity 



li^go 
-1 



{goEh - /I' 
1 



flbgi?* 1 — Cbg/c 



(67) 
(68) 



1 



1 Airfi^/m 

a go - A'^/Eb by^ 



which can be positive or negative. The ratio go/ A is also 
easily eliminated with the identity 



(59) 



AttA^ 



— R*a,i, 



(69) 



The large Ei, limit of this expression is the inverse of the 
background scattering length: 



Our final expression for the scattering amplitude is thus 



1 47rft2/m 

flbg go 



1 



60f 



(60) 



1 

Jiq 



2.2 
,-q b 



1 - (1 - a/flbg) 



The location of the Feshbach resonance corresponds to 
an energy Ej^ of the closed-channel molecule such that 
1/a = 0, so that 



I.e. 



1 

fk 



^k^b-' 



1 - (1 - a/flbg) ■ 



Et = 



A^ 



go + 4:7T^/'^h'^b/m' 



(61) 



Supposing that the energy Ei, of the closed-channel molecule 
varies linearly with the magnetic field B, 



We deduce the effective range in terms of a, Obg, R* and b 
Te = —2R^,(l — abg/a) H — -1= . 



- q erf (qb) — q 
(70) 

-ifcerf {^ikb)+ik. 

(71) 
id 5: 

(72) 



EbiB) ^ E^ + ^itiB - Bo), 



(62) A. 2 Conditions for reaching the zero-range limit 



and introducing the (possibly negative) resonance width 
AB such that 

/12 

fibAB El (63) 

.90 



In a degenerate gas, we expect that the zero-range limit is 
reached and the equation of state E{l/a) becomes model 
independent provided the conditions 



we get from (|59l) exactly the B-field dependence PH)) for 
the scattering length a in free space. 



kb^ 1 
1 _ 1 
fk a 



ik 



(73) 
(74) 
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hold for k ^ kp and for k ~ 1/a on the BEC side [38] , 

Wc now assume that (|73|) holds, and we derive validity 
conditions for (|74p in the case of our two-channel model. 
Following the usual procedure, we use the second order 
expansion 



1 1 ■, 1 

= — \- tk Tp 

fk a 2" 



(75) 



which reduces to ((74|) if the effective range is small enough 
to have 



2' 



\e<. 



-+ik 
a 



(76) 



Taking into account the condition ([75)) as well as the ex- 
pression (|72[) of Te for the two-channel model, and using 
\i + l/{ka)\ > 1 and |1 -I- ika\ > 1, we find that the condi- 
tion ([76ll becomes equivalent to 



k'^R, 



(,_^y 



- +ik 
a 



<1. 



(77) 



A second condition is obtained by imposing that the higher 
order terms ... in (j75[) are indeed negligible compared to 



;|/2. For this to hold we need not only ([75)) . but also 



k' 



fc^«|Q^|i.e. 



fc^|abg|i?* 



1- 



abg 



<1. 



(78) 



We conclude that to reach the universal zero-range 
limit, it is sufficient to fulfill the three conditions (|73l77l78p . 
All three conditions are satisfied for the experimental data 
on ^Li from [5j reproduced in FiglU where l/kp > 250 nm 
and \a\ > 200 nm: using the values |abg| = 74 nm and 
i?* = 0.027 nm (see section [2731) ^^.d taking the estimate 
6 ~ 2.1 nm obtained by imposing that the effective range 
on resonance ([65]) should be equal to the value 4.7 nm 
obtained in the multi-channel calculation of [87^ , we find 
that the left hand sides of Eas. (|73l77l78p are respectively 
smaller than 10-^ 3 • 10"'* and 7 • 10"^ 

To be complete, let us point out that the conditions 
(J77I78P are not stricly necessary since it is possible that 
the ... in ([75)) are not neghgible compared to k^\re\/2 
and ([71)1 nevertheless holds. We have chosen these condi- 
tions because they simplify the discussion, in particular 
they automatically ensure that the analytic continuation 
fig has a pole at g ~ 1/a so that the dimer energy is 
~ —lT?/{ma?). For this last property, in the regime b -^ a, 
one directly gets from ()70p the truly necessary condition 



|(l-a/abg) ^ + (l-abg/a) ^a/i?J > 1. 



(79) 



Close to the Feshbach resonance, a ^ |abg|, this reduces 
to the simpler condition 

a:$>R^. (80) 



B Effect of a difference between the trapping 
potentials seen by a closed-channel molecule 
and by an open-channel pair of atoms 

In this appendix, we consider the physical consequences of 
the fact that, in general, the trapping potential Ub{r) seen 



by a closed-channel molecule is not simply the trapping 
potential 2 C/(r) seen by two separated atoms in the open 
channel: 

Ub{r)^2U{r), (81) 

and we derive conditions for the gas to behave approxi- 
mately as if we had Ub{r) = 2U{r). Note that a signifi- 
cant deviation in ([5T)) is not a priori excluded for the laser 
traps used e.g. in [51|56] for lithium, with a laser wave- 
length Al ~ 1 /Ltm: Contrarily to the atomic lightshift, the 
lightshift of the closed-channel molecule may be sensitive 
to the laser detuning from transitions to bound states of 
the 1^E^{A) Born-Oppcnheimcr interaction potential be- 
tween a ground-state atom and an excited atom, consid- 
ering the fact that the minimum of this excited-ground 
interaction potential is separated from the dissociation 
limit of the ground-ground interaction potential by only 
~ 0.6 ^m^^ < 1/Xl in spectroscopic units [88| . 



B.l Position dependent scattering length and effective 
magnetic field 

At the two-body level, ([5T)) leads to an effective posi- 
tion dependent scattering length for two opposite spin 
atoms. To define such a local scattering length around 
point ro, we consider the "tangential" free space problem 
where the trapping potentials are replaced by the con- 
stant values Ub{ro) and U{ro) and we call a(ro) the re- 
sulting scattering length. We simply have to revisit the 
calculations of appendix |^ In ([5^ one has to replace 
Et by Eb + Ub{ro) to include the trapping energy shift 
of the closed-channel molecule; similarly in ([55)1 one has 

to replace ^^ by ^^ -I- 2[/(ro) to include the trapping 
energy shift of the two atoms; finally, for atoms with in- 
coming wavevectors ±ko, the scattering state eigenenergy 

is now E = — ^ + 2U{rQ). The net effect on the scattering 
amplitude fkg in (|56|) is to apply the substitution 



Eb 



Eb + Ubiro)-2Uiro), 



(82) 



so that a(ro) is simply obtained by applying ([5^ to 

A powerful interpretation of this substitution proce- 
dure may be given for a magnetic field independent fib'- 
We can rewrite the term Eb{B) + Ub{r) appearing in ([T]) 

as £'6[i?cff (r)] + 2 U{r) provided we set 



Bcff(r)=S 



Ub{r)-2U{r) 



Mfc 



(83) 



Thus we can replace Ub{r) by 2U{r) in (|T]) provided we 
also replace the homogeneous magnetic field S in (|T]) by 
the inhomogeneous effective magnetic field Bcff{r). One 
then immediately obtains the position dependence of the 
scattering length. 



a(r) = afs[Bofr(r)], 



(84) 



the function afs(i3) being given in ([lOp . 

This effective magnetic field picture is also very efff- 
cient at the many-body level. In the grand canonical point 
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of view, one adds the term —fiN to the Hamiltonian _ff2, 
where /i is the chemical potential, and the particle number 
operator N is defined as iV = Nf + 2Nb, Nf giving the 
number of fermions and Nf, the number of bosons, since 
this is the quantity which commvites with H2 ■ In the local 
density approximation, in the general case Ub 7^ 2C/, we 
thus see that the gas behaves locally like a free space, ho- 
mogeneous gas of chemical potential /i — t/(r), in presence 
of the locally homogeneous magnetic field Bas{r). 



B.2 Shift of the resonance location 

We now consider the case of an optical trap such as the 
one at Rice [5]: U{r) and Ub{r) are the lightshift potentials 
experienced by an atom and by a closed-channel molecule. 
Being both proportional to the laser intensity, they are 
mutually proportional: 



Ub{r)=2{l + f^)U{r) 



(85) 



where the dimensionless parameter 77 is a convenient mea- 
sure of the relative deviation between Ub and 2U. At large 
r the laser intensity vanishes, so do Ub and U. 

In practice, the gas is trapped around the intensity 
maximum of the laser field, located in r = 0, and the ex- 
tension of the cloud along the axial (resp. radial) direction 
is small compared to the Rayleigh length (resp. the waist) 
of the laser beam; thus, to leading order, U{r) ~ U{0), we 
can temporarily forget about the spatial variation of Beg 
and take BcS ~ _B — SBq, where 



5Bo 



-Ub{0) + 2U{0) 

^J.b 



2riU{0) 



^Jib 



(86) 



Since the free space scattering length afs(-B) only depends 
on the difference B — Bq, replacing B by Bcs is equiva- 
lent to shifting the position _Bo of the Feshbach resonance 
by the quantity 5Bo. Introducing the length _R* defined 
in (im, the ratio of this shift to the width AB of the 
Feshbach resonance may be written as 



5Bo 
AB 



2mU{0) 
fi2 



flbgi?*?]. 



(87) 



In the Rice experiment, the trap depth — C/(0) is at most 
ks X 2.71 ^K [52], so that we have 



\SBo\ 
\AB\ 



< 



7500 



where we have used the values of Obg and -R* given below 
([^^ and the mass of lithium. Except if 77 takes truly large 
values for the laser wavelength Al = 1064 nm used in [S], 
we expect only a small shift of the resonance location due 
to the laser trap, relative to the resonance width. 

In [5] the value of the scattering length as a function of 
B was calculated essentially from the free space formula 
a[s{B) ignoring the shift SBq due to the trap, and taking 
for the resonance location Bq the value measured in [SS] in 
a laser trap with different laser wavelength Al = 1030 nm 



and intensity. On the contrary, our definition of a as a 
function of B includes the shift SBq, a = afs(i? — SBq), as 
already defined in ([9]). This may introduce a systematic 
discrepancy between our theory curves and the experi- 
mental data in Fig[l] in the form of shifts of the abscissas 
of the experimental points. We thus define the horizontal 
shift 



1 

kpa 



1 



1 



kpUfsiB - SBq) kpafsiB) 
SBq 1 (1 - flbg/a)' 



AB kpabg 1 - (1 



ojSBq/AB 



(89) 
(90) 



To leading order in rj (that is for \SBq/AB\ <C 1), and 
using the fact that the data from [5] reproduced in Fig. [1] 
is sufficiently close to resonance to have |1 — abg/a| < 1-4, 
this simplifies to: 

6 (A\ ^ ^ 1(1 _ «,,/„)2. (91) 

ykpa) AB kpay^g^ "' ' ^ ' 

Using the parameters of the Rice experiment [89^ and 
Eq.dHlI) this gives 



1 
kpa 



< 



1000 



(92) 



Even without knowing 77, we can thus resonably expect 
that this shift is negligible on the scale of Fig. [2 i.e. that 



1 
kpa 



<1. 



(93) 



This is also consistent with the fact that no systematic 
horizontal shift is apparent in FiglT] between the experi- 
mental data and the theory curves. 



B.3 Effect of the position dependence of Bcs{r) 

We now discuss the effect on the many-body properties 
of the gas of the inhomogeneity of Bcs{r), that is of a 
spatial dependence of the scattering length a(r), see ([M|) . 
now taking into account the r dependence of U{r). We 
rewrite (1831) as 



Scff(r) = B- 



2TjU{r) 



Aifc 



= B-SBq- 



27?[[/(r) - U{Q) 

^J■b 



(94) 



When inserted in ((84)) this leads to the following position 
dependence for the scattering length: 



1 



1 



a(r) — Cbg a — flbg 



77^[[/(r)-f/(0)], (95) 



where we recall that a is the scattering length in the trap 
center. Right on resonance, \a\ = -I- 00, i?* is one of the 
contributions to the interaction effective range, see (|65|) . 
so that the zero-range limit implies fc/^li?, | <C 1. Since the 
typical value of C/(r) — t/(0) over the cloud size is of the 
order of H'^kp/2m, we see that the small parameter kpR^, 
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appears in (|95p so that we may expect that the position 
dependence of a(r) is negligible in the zero-range hmit (for 
a fixed value of 77) . 

Let us investigate this expectation more quantitatively 
and beyond the unitary limit, however still in a limiting 
case with simplifying assumptions. First, we assume that 
we are sufficiently close to resonance to have |a(r)| 3> |abg| 
over the atomic cloud size. This implies that the excur- 
sion of BcS (r) over the cloud size is small as compared to 
|zii?|; since the excursion of U{r) over the laser trap size 
is at most |C/(0)|, we note that this condition is automat- 
ically satisfied for \SBo\ <C \^B\. Second, we shall assume 
that the gas may be treated in the zero temperature lo- 
cal density approximation. Third, we use the quadratic 
approximation (J13p for U{r), further supposing for sim- 
plicity that the harmonic trap is isotropic, a case to which 
one can always reduce within the local density approxi- 
mation by a rescaling of the coordinates, see appendix [Cl 
Equation ([55)1 then reduces to 



a{r) 



1 

--ry- 

a 






2 2 
muj r , 



and the gas density profile n(r) is given by 



li-U{0) 



-mu?r^ + i^°'^ 



71 (r). 



<^y 



(96) 



(97) 



where /i''°'°[7i, 1/a] is the chemical potential of a free space 
homogeneous gas of density n and scattering length a. 

Let us derive simple conditions to have a weak effect of 
the scattering length inhomogeneity on the density profile. 
Some of these conditions will be expressed in terms of kp 
defined by (fU]) . 

In the BCS regime, a < and kp\a\ < 1, we have 
fi - U{0) ~ ksTp = h^k\/2m. Over the Thomas-Fermi 
radius R of the cloud, such that ^muj^R^ ~ ksTp, we see 
that the relative change of l/a(r) is negligible if 



|1 - a/a{R)\ ~ {kF\a\){\r^\kFR*) < 1. 



(98) 



This is satisfied in the zero-range limit, for fixed 77, accord- 
ing to the condition ([77| written for k — kp. 

In the unitary limit kp\a\ 3> 1, the previous reasoning 
fails, because requiring that the relative change of a{r) is 
small is too stringent. Taking for simplicity the limiting 
case \a\ = -f 00, it is sufficient to check that the gas remains 
locally unitary except in a small region near the edge of 
the cloud. First, to zeroth order in 77, we calculate the 
local Fermi wavevector fc^°™ (r) using the unitary equation 
of state fi^°'^[n, 1/a = 0] = ^n^{k^°'^f /2m, where k^""" 
is defined in (fT7|) . and using the value of the chemical 
potential /i — C/(0) = ^/^ksTp. We obtain 



jLhom 



^') = 7rrA^-m 



^1/4 



i?2 



1/2 



(99) 



where R is the unperturbed Thomas-Fermi radius such 
that muj'^R^ /2 = y^ksTp. Then, we include the 77 posi- 
tion dependent term in l/a{r) to calculate at which dis- 
tance I from the edge of the cloud the gas starts leaving 



the unitary regime, that is 

k'P"^iR-l)\a{R-l)\ = l. (100) 

Assuming a priori I ^ R we obtain from (|99p 
I ^3/2 



R 



-ivkFR*r 



(101) 



We conclude a posteriori that l/R ^ 1 and that the spatial 
dependence of 0(7^) is negligible when 



|77|fcFi?* < 1. 



(102) 



As in the BCS case, this is equivalent for a fixed I77I to the 
zero range condition (|77p written here with k ~ kp and 
\a\ — -l-oo. 

Finally we turn to the BEC regime a > and kp\a\ < 
1. We then use the approximate equation of state 



fi^°^^[n{r),l/a{r)] 



-Kh' 



2ma?'{r) 2m 



aa{r)n{r) (103) 



where a ~ 0.60 is the proportionality factor in (|20p be- 
tween the dimer-dimer scattering length ad and a. In par- 
ticular, this requires a :^ R^, see ([50]) . Replacing l/a(r) 
by ([M]) in the above expression, we find that the leading 
term in l/a'^{r) gives rise to an external potential that 
sums up with the trapping potential to give 



U,s{r) = U{r) - 



2ma'^{r) 



(104) 



= UMO) + i-,2ffr2 - ^^^^(mc.V2)2(105) 

We see that the atomic oscillation frequency is renormal- 
ized into 



^cff 



1 



2riR, 



1/2 



(106) 



To have a small effect of the scattering length inhomo- 
geneity on the trapped gas, we first require tOcS ~ cj so 
that 

|77|i?Ja<l, (107) 

a condition with no counterpart in the BCS regime and 
more stringent than (|102p since kpa < 1 here. Again, we 
find that we recover, for a fixed I77I, the zero range condi- 
tion (|77p written here for k = 1/a. Second we require that 
the r^ term in Ucs{t) remains small as compared to the 
r^ one over the Thomas- Fermi radius R of the condensate 
of dimers, calculated here to zeroth order in 77; we reach 
the second condition 

!I^^n^^^2R^ ^ {^kpR^f U-^kpaX ' « 1. (108) 

Since kpa < 1 here, this second condition is automati- 
cally ensured by the first condition (jl07p . Third, we re- 
quire that the relative variation of l/a(r) is small over the 
Thomas-Fermi radius R of the condensate, to ensure that 
the spatially inhomogeneity of a(r) in the mean field term 
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[irh? / (2'm)]aa{r)n(r) has a small impact on the density 
profile: 

/ 5 \ 2/5 

\l-a/aiR)\ ~ {kFa){\r,\kFR*) ( ^^f«d j « 1. (109) 



which is also implied by (|107p . 

To summarize, discussing perturbatively the effect of 
the position dependence of Bcs{r) on the gas density pro- 
file in the local density approximation, we find close to 
the Feshbach resonance (|a| ^ kbgl, and also a 3> -R* in 
the BEC regime) that the effect is small under the con- 
dition ([98| in the BCS regime, under the condition ()f02|l 
in the unitary regime, and under the condition (|f07p in 
the BEC regime: In all three cases, this is the zero range 
condition (|77p . written respectively for k — kp,k = kp 
and k — I/a, and multiplied by the parameter \ri\. If \r]\ is 
not too large, neglecting the r dependence of the scatter- 
ing length is thus automatically justified in the zero-range 
limit. 

Interestingly, in the unitary limit, it is possible to re- 
late rjkpR* to the quantity 6{l/kpa) introduced in the 
previous subsection: 



IvlkpR*^ 



1/2 



kpa 



U{0) 



(110) 



Experimentally, the gas is confined at the bottom of the 
laser trap so that [n — U{0)]/\U{0)\ <^ I (its value near 
unitarity at Rice is 0.15 [89j). Together with the assump- 
tion ([M]) this implies that the right-hand-side of (jIIOp is 
^j; I, and that the position dependence of a(r) has a neg- 
ligible effect on the gas density. 

A quantitative extension of the present qualitative dis- 
cussion is to evaluate the density change Sn{r) of the gas 
to first order in 7y. To this end, one simply expands (j97p 
to first order in 77 around the solution no{r) with 77 = 0, 
writing to first order n{r) = nQ(r) + 6n(r) and fj, = fj,o + Sfi: 



dn{r) 






(111) 



where we have introduced the variation of I/a(r) around 
77 = to first order in 77, not supposing \a{r)\ 3> jobgl in 
13: 



<5 ( ^ ) (r) 



I 



ab: 



rj — T^rriLd r . 



(112) 



(5/i is then determined by the condition of a fixed number 
of particles /^,p (PrSn{r) = 0, where Rq is the Thomas- 
Fermi radius of the gas for 77 = 0. In this way, taking 
the variation of p30p to first order in 77 around 7; = 0, 
one can calculate the effect of a non- uniformity of a(r) on 
dE/d{l/a) to first order in 77. We give here the result for 
the unitary gas: 



6 [dE/d{l/a)] 
dE/d{l/a)\r,=o 
3157r 
"''2048 



= rjkpR^ 



/_6£ _637r\ 1/4 
^257r o'^R 1^ '' 



256 



■e^'k 



Fa.be 



(113) 



where we recall that ^ = /(O), C = -/'(O) , see dH]), and 
/*^2)(o) is the second order derivative of f{x) with respect 
to a; in a; = 0. Note that, in this first order variation with 
respect to 77, the trapping potential U{r) is kept fixed so 
that kp is also fixed. 

The calculation (|113p shows that, at unitarity, a mea- 
surement of the variation of the number of closed-channel 
molecules Nb as a function of Ub{r)/U{r) may give access 
to Z'^^) (0). It also shows that, in the unitary limit, the con- 
dition to neglect abg with respect to a(r) as done in the 
qualitative reasoning around (|I02p . is fcF|abg| <C 1. 



B.4 Metastability issues for a position dependent 
scattering length 

We now discuss whether it is energetically possible for a 
pair of atoms to form a dimer which could escape from the 
trap by tunneling. The energy variation in such a process 
is 



AE 



mo?' {00) 



-2^ 



(114) 



where /x is the chemical potential of the trapped gas and 
a{oo) is the effective scattering length at infinity, which we 
assume to be positive. We discuss under which conditions 
this process is energetically forbidden, i.e. AE > 0. 

We assume that the gas has the usual Thomas-Fermi 
density profile and is weakly affected by the position de- 
pendence of i?cff (r), see subection lB.31 We distinguish be- 
tween two regimes for the scattering length a in the trap 
center: the negative-a regime where 1/a < 0, and the BEC 
regime where < kpa < 1. Under the usual assumption 
that the gas is confined in a small region around the bot- 
tom of the laser trap we have fi ~ U{0) in the negative-a 
regime and /i ~ U{0) — ti^/{2ma'^) in the BEC regime. 
We also make the simplifying assumption a{oo) ^ |abg| 
so that (l95]l reduces to 



1 1 in 

- + 2r]—RMiO). 



a(oo) a h? 

In the BEC regime one easily deduces 

AE 



2\um 



l + A-B 



where 



„ R* 
A= 2ri^ 
a 

B^2rj'^\Um; 



(115) 
(116) 

(117) 
(118) 



since we already assumed (J107I) we automatically have 
\A\ <C 1. In the negative-a regime, Eq. (|115p immediately 
gives an upper bound on l/a{oo), from which we get 



AE 

2\um 



>l-B. 



(119) 
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In both regimes we conclude that under the condition The variation of the gas energy is then 



\B\ « 1 



(120) 



we have AE > 2\U{0)\ so that the process is energetically 
forbidden. 

Finally we note that under the assumption |a| S> |abg| 
we have from (|90l87p 



S[ - — - vkpR* ■ ^ 



so that 



B'- 



kpa 






(121) 



(122) 



In a typical experiment such as the one at Rice we have 
ksTp < \U{0)\; thus the condition \B\ ^ 1 is satisfied as 
soon as \S{l/kpa)\ <C 1, which is likely to hold at Rice, 
see (I92l). 



C Calculation of dE/d{l/a) in the local 
density approximation 

Restricting to a spin balanced gas, we first consider the 
general case of the two-channel model ([T]), in which Ub{r) ^ 
2C/(r). As explained in the appendix [Bl where an effective 
position dependent magnetic field _Boff(r) is introduced, 
this leads to a position dependent scattering length a(r). 
As a consequence, in the local density approximation, here 
at zero temperature, the gas density n(r) is given by 



^ = C/(r)-|-^'^°-Kr),lMr)] 



(123) 



where /J,^°™[n, 1/a] is the chemical potential of the free 
space homogeneous gas of density n and scattering length 
a. The mean energy and particle numbers in the gas are 
given in the same approximation by 

E= /"dV{C/(r) + e''°"^[n(r),l/a(r)]}n(r) (124) 
N = I £rn{r) (125) 

where e'^°™[n, 1/a] is the mean energy per particle of the 
free space gas. We wish to calculate the derivative of E 
with respect to 1/a, where a is the scattering length in 
the trap center r = 0, with the constraint that N should 
be fixed. 

Let us perform an infinitesimal variation 6 (i) of 1/a. 
This leads to a variation Sfi of the chemical potential and 
6n{r) of the density, that may be calculated from (|123p but 
this is not required here. This also leads to a variation of 
l/a(r), according to ([95]) : 



1 



a(r) 



1 - Obg/Q(r) 
1 - flbg/a 



(126) 



6E = / d-'r Sn{r) 



-n{r) 



U{r) 



^hom I 



7i(r),l/a(r)] 



de 



hoin 



dn 



[n(r),l/a(r)] 



a(r) 



n^hom 

n{r)^r^—Jn{r),l/a{r)]. (127) 
d{l/a) 



Writing ehom[n, 1/a] = ^''°'"[iV, V, l/a]/N, where V is the 
volume and N the particle number of the homogeneous 
gas, we find the relation 



g^ho 



dn 



n, 1/a] = n''°"'[n, 1/a] - e''°'"[n, 1/a]. (128) 



Using this expression and (|123p one finds that a simpli- 
fication occurs in (|127p . Furthermore, the fact that the 
particle number is fixed imposes 



(fr5n(r) = 



(129) 



so that (|127|) leads to 



dE 
d{l/a) 



N 



= / d'^rn{r) 



1 - abg/a(r) 



g^ho 



d{l/a) 



flbg/a 

[n(r),l/a(r)]. (130) 



We note that an alternative derivation of this result is 
possible: One can introduce the density of closed-channel 
molecules njj°'"[n, 1/a] for a free space gas of density n 
and scattering length a, a quantity that can be related to 
the derivative of e'^°™[n, 1/a] with respect to 1/a simply 
by applying the general formula ([8]) to a free space system. 
Then setting for the trapped system in the local density 
approximation Nt = J d'^rn^°"^[n{r),l/a{r)] and using 
([8]) for the trapped system leads to (|130p . 

In practice, in the present paper, we perform explicit 
calculations in the limiting cases presented in W2.'2[ One 
takes the zero-range limit, which also allows to neglect 
the spatial variation of a(r) under conditions discussed 
in the appendix [BJ Then e^°'^ is given by the relation 
()16|) involving some universal function f{x). For the free 
space chemical potential we then introduce the convenient 
parametrization 



fJ- 



hom 



71^ 1/a] 



2ma^ 



Uhom^ 



(131) 



where fcj^°™ — (37r^n)^/'^ is the ideal gas Fermi wavector, 
and the function u is 



u(x) 



fix) - Ixf'ix) 



(132) 



One also restricts to a harmonic trap for C/(r), see ((Tc 
The surface n(r) = is then an ellipsoid with Thomas- 
Fermi radii Ra along the directions a. 
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For a < we set 



2„2 



fJ- 



2m 



(133) 



with q > 0. Then uJaRa — fiq/m, and the density depends 
only on X > such that X'^ = J^a ^a/^Q- Equation p^ 
then reduces to the impUcit equation 



q^a'^{l-X'^) = u 



1 



k!f"^{X)a 



(134) 



that one can solve numerically, using the interpolation for 

f{x) given in appendix ID] 

For a > 0, /^'^°'"[n, 1/a] tends to -h^ / {2ma'^) forn -> 

so that we set 

h^q^ fe , , 

A* = -TT- - T. 7 135 

^ 2m 2ma2 ^ ' 



with g > 0. Then again LUaRa ^ hq/ni, and (|123p reduces 
to 

1 



q^a\l-X^) = 1 + ' 



k'f'^{X)a 



(136) 



to be solved numerically in general. 

One can then calculate dE/d{l/a) and kp defined in 
(fT4|) from the expressions valid on both sides of the reso- 
nance: 



dE 3Nh^IodXX^k^°^{Xrf 



°(X)a 



d{l/a) 
kp 



IOto 

7ri/3 



An 



J^ dX X^kf'^{X)3 

1/6 



fl37) 



dXX^kf"" 



{Xf 



(138) 



The analytical calculations performed in the weakly inter- 
acting regimes are also obtained from the above formulas. 



D Interpolation formula for the homogeneous 
gas energy 

We explain how an interpolation formula was constructed 
for the function / defined in (flB]) , which was then used to 
produce the solid line in Fig[T^ and FigdJD. 

On the BEG side of the resonance, we applied a cu- 
bic spline interpolation to the points in Table I of [32] 
with values of x = l/(fc^°™a) equal to 0, 1, 2, 4 and 6, 
using the last column of Table I (the one said to give 
E/N — eh/2) and adding to the corresponding data the 
quantity —fi?/{2ma?). We added a home made point at 
X = 8 where the value of f{x) and its first order deriva- 
tive /'(x) are obtained from the bosonic Lee-Huang- Yang 
expansion p^ . 

On the BCS side of the resonance we used the empirical 
interpolation formula of [90] . 

fs{x) ~ ai — 0:2 arctan[a3x(/3i — x)/((32 — x)]. (139) 

Constraints on 01,0:2,03 and /3i,/?2 are obtained from 
the values /(O), /'(O) and from the large |x| expansion 



f{x) = \ + ci/x + C2/x^ + 0{l/x^), where the coefficients 
ci and C2 are given in ([22]) . This differs from [QO' where the 
coefficient C2 was not used and /32 was deduced from a fit to 
the Monte Carlo data of [32] on the BCS side. Another dif- 
ference with [5D] is that we took the value /'(O) = —0.95, 
that we extracted from [34] using the relation ((34)) . Also, 
we set the value of /(O) to 0.41, which is within the error 
bars of [32j . in order to reduce the spurious discontinuity 
of the second order derivative of f{x) in .t = between the 
interpolation formula for x < and the spline for x > 0. 
Ah this leads to oi ~ 0.41, 02 ~ 0.3756, 03 ~ 1.062, /3i ~ 
0.9043, /32=i 0.3797. 

A test of this interpolation formulation is to compare 
its prediction for /'(—I) with the value ~ —0.13 that one 
can extract from [34] using ((34]) ; one gets the satisfactory 
result /s(-l) ^ -0.14. 
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